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Random matrix theory (RMT) is a powerful statistical tool to model spectral fluctuations. In addition, RMT 
provides efficient means to separate different scales in spectra. Recently RMT has found application in quantum 
chromodynamics (QCD). In mesoscopic physics, the Thouless energy sets the universal scale for which RMT 
applies. We try to identify the equivalent of a Thouless energy in complete spectra of the QCD Dirac operator 
with staggered fermions and SUc{2) lattice gauge fields. Comparing lattice data with RMT predictions we find 
deviations which allow us to give an estimate for this scale. 
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In recent years, RMT has been successfully 
introduced into the study of certain aspects of 
quantum chromodynamics (QCD). The interest 
focuses on the spectral properties of the Eu- 
clidean Dirac operator. For the massless Dirac 
operator _^[f7] with staggered fermions and gauge 
fields U S SUc{2) we solve numerically for each 
configuration the eigenvalue equation 
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The distribution of the gauge fields is given by 
the Euclidean partition function. Examples of 
the spectra are shown in Fig. |l|, where the av- 
erage level densities for 16'* and lO"* lattices are 
shown. It should be pointed out that we have 
V/2 = 32768 resp. V/2 = 5000 distinct positive 
eigenvalues of each configuration, so that there 
are millions of eigenvalues at our disposal. 

As the gauge fields vary over the ensemble of 
configurations, the eigenvalues fiuctuate about 
their mean values. Chiral random matrix the- 
ory models the fluctuations of the eigenvalues in 
the microscopic limit, i.e. near A = |0] as well 
as in the bulk of the spectrum Our main 
question is to what scales RMT does apply in 
QCD. In disordered systems the Thouless energy 
Ec determines the scale in which fluctuations are 
predicted by RMT. Beyond this scale deviations 



0.2 

P(A) 
V 

0.1 



1 1 1 1 1 1 1 
" (3 = 2.4. L = 


1 1 1 1 1 1 1 1 1 1 
16 ^ ^ 


1 1 


921 conf. 


'^y^ — ^ 






^ \ 






1416 conf. 1 




^ 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 


J 1 



*Poster presented by M.E. Berbenni 



1 2 3 4 5 
X 

Figure 1 . Spectral density of the lattice Dirac op- 
erator for A > for two different lattice volumes 
and gauge coupling values. 

occur. In QCD an equivalent of the Thouless en- 
ergy is Armt In the microscopic region it 
scales as 

Armt/^? (X W, (2) 

V is the lattice volume, D is the mean level spac- 
ing. As argued in jj] a corresponding effect should 
also be seen in the bulk of the spectrum. 
The staircase function N{X) gives the number of 
levels with energy < A. In many cases it can be 
separated into 

A^(A) = A^avc(A)+iVfluc(A). (3) 
A^avc(A) is determined by gross features of the 
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Figure 2. Number variance (left) and spectral rigidity (right) on a 16'* lattice. The solid lines are RMT 
predictions, the dots numerical data. On this scale the data points do not depend on unfolding. 



system. Naud^) contains the correlations to be 
analyzed. RMT makes predictions for the fluc- 
tuations on the scale of the mean level spacing. 
The influence of the overall level density must 
be removed by numerically unfolding the spec- 
tra through the mapping Xi — s- Xi = iVavc(Ai). 
For the new sequence we then have Nnvoix) = x 
i.e. the mean level spacing is unity everywere 
1/Pavc(a;) = 1 where Pavo(a;) = dN^yc{x)/dx. The 
extraction of 7Vavo(A) is highly non-trivial, be- 
cause little is known analytically about the level 
density of QCD spectra. 

However, there are several phenomenological 
unfolding procedures, e.g. ensemble unfolding, 
where one divides the energy range in m bins 
of width AA and averages the density p(A, A -|- 
AA) for each bin over all configurations. Then 
the staircase function Afave(A) = p{^i:^i + 

AA)AA, with Am = A is calculated. Furthermore 
there is configuration unfolding, where N{X) is 
fitted for each configuration to a polynomial of 
degree n. Strong couphng expansions for SUc{2) 
with staggered fermions [pi and 1/Nc expansion 
of the QCD level density motivate this ansatz. 
For technical details and further unfolding proce- 
dures see 0. Whatever approach one uses, the 
mean number of rescaled levels in an interval of 
length L in units of the mean level spacing should 
equal L. This assures that the unfolded spectrum 
has mean level density unity. 

We compared RMT predictions for two-point 
correlators with lattice data for two quantities. 
First, the level number variance, which measures 
the deviation of the number of eigenvalues na{L) 



in an interval [a, a + i] from the expected mean 
number L 

Y?{L)^ < {L-n^{L)f >. (4) 

< • • • > is the spectral average, (• • •) the ensem- 
ble average. Thus, an interval of length L con- 
tains on average L± y/iy^L) levels. For uncorre- 
lated Poisson spectra S^(L) = L. RMT predicts 
stronger correlations: S^(L) ~ \ogL. The second 
two-point correlator we considered is the spectral 
rigidity, defined as the least square deviation of 
iV(A) from the straight line 

ML) = <^|min^,B^deV(e) -A^- Bf"^. (5) 

For this quantity RMT predicts A3(L) ~ logL. 
In Fig. 1^ the RMT results for these statistical 
measures are compared with lattice data. The 
wealth of data allows us to analyze higher order 
correlations. Again we see good agreement (see 
Figs. 11 and 12 in 0). 

With ensemble unfolding of the data we obtain 
for S^(L) the curves plotted in Fig. |[ Indepen- 
dently of the spectral region considered and of 
(3 we find that the point where the deviation sets 
in, scales as 

Xkmt/D ~ O.SVV. (6) 

This should be compared with the result obtained 
in for the microscopic region 

Armt/^^0.3...0.7VV. (7) 

With polynomial unfolding the scaling law (^ 
vanishes. The deviation point appears to be the 
same for different lattice sizes (Fig. ^) . 
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Figure 3. Deviations from RMT predictions for different lattices sizes V and gauge couplings /3. Shown 
are different regions of the spectrum as indicated in the upper left part of the plots. 
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Figure 4. Deviations from RMT predictions with 
configuration unfolding of the data for different 
lattice sizes and /3 values. 

In order to find out, if these deviations are 
due to a Thouless energy, we performed a Fourier 
analysis of the oscillations of the staircase func- 
tion . From it we concluded that the deviations 
of the data obtained with polynomial unfolding 
are due to a non-polynomial-like part in the aver- 
age level density and not to an equivalent of the 
Thouless energy. 

In conclusion, analyzing some of the statistical 
properties of complete eigenvalue spectra of the 
Dirac operator for staggered fermions and SUc{2) 
gauge fields for various couplings and lattice vol- 
umes, we find the scaling behavior of the equiva- 



lent of the Thouless energy. Using ensemble un- 
folding, we have Xrmt/D = CVV. The con- 
stant is approximately C « 0.3 which is compat- 
ible with the result obtained in |^ for the micro- 
scopic region of the spectrum, where the scaling 
(|^) was found. By unfolding each configuration 
separately, we do not see any scaling of this type. 
Hence the Thouless energy is due to fluctuations 
in the ensemble. 
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